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ABSTRACT 

As part of a survey for radio pulsars with the Parkes 64-m telescope we have discovered PSRJl 622-4950, a 
pulsar with a 4.3-s rotation period. Follow-up observations show that the pulsar has the highest inferred surface 
magnetic field of the known radio pulsars (B~3 x 10^"* G), exhibits significant timing noise and appears to have 
an inverted spectrum. Unlike the vast majority of the known pulsar population, PSRJl 622-4950 appears to 
switch off for many hundreds of days and even in its on-state exhibits extreme variability in its flux density. 
Furthermore, the integrated pulse profile changes shape with epoch. All of these properties are remarkably 
similar to the only two magnetars previously known to emit radio pulsations. The position of PSRJl 622- 
4950 is coincident with an X-ray source that, unlike the other radio pulsating magnetars, was found to be in 
quiescence. We conclude that our newly discovered pulsar is a magnetar - the first to be discovered via its radio 
emission. 

Subject headings: pulsars: individual (IE 1547-5408, PSR J1622-4950, XTE J1810-197) — stars: magnetars 
— stars: neutron 



1. INTRODUCTION 

Magnetars are slowly rotating neutron stars that, in contrast 
to ordinary pulsars, are not powered by their spin-down en- 
ergy losses, but by the energy stored i n their extremely large 
magnetic fields, typically > 10^"^ G (iDuncan & ThompsonI 
[1992). They are divided into two groups. Soft Gamma-ray 
Repeaters (SGRs) and Anomalous X-Ray Pulsars (AXPs), 
whose relation to each other and to other kinds of neutron 
stars is stiU a matter of debate (iMereg hettill2008h . SGRs 
were discovered as sources emitting short repeating bursts 
in the hard X-ray/soft 7-ray band, while AXPs were first 
detected in the soft X-ray band. Recently, sources have 
been found that show properties of both groups, suggesting 
that AXPs and SGRs belong to the same family ( React al. 
120091: iMereghetti et al.l l2009h . Many attempts have been 
made to find pulsed r a dio emission from magnetars (see 
e,g^_^ lKrisseran 119851: iGaensler et all 120011: iBurgav et al.l 
I2006h but it was not until 2006 that the first detection 
was reported (Camilo et al. 2006). To date only two AXPs 
have confirm ed radio pulsations, namely X TE Jl 8 10-19 7 
(ICamilo et al.ll2006.) and IE 1547-5408 (Camil oetani2007ah . 
Both of these sources are transient AXPs, a subgroup of mag- 
netars that occasionally undergo very bright X-ray outbursts. 
Their pulsed radio emission was discovered in connection 
with these outbursts, and is observed to fade with time after 
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an outburst. 

Here we report on the discovery of PSR J 1622-4950, the 
first magnetar discovered in the radio band, and discuss its 
properties in relation to the two magnetars previously known 
to emit radio pulsations. 

2. OBSERVATIONS AND HISTORICAL DATA 

2.1. Parkes Observations 

The High Time Resolution Universe (HTRU) survey for 
pulsars and fast transients is cu rrently being c arried out at 
the Parkes 64-m radio telesc ope (iKe ith et a l.ll201Q.), usin g the 
20-cm multibeam receiver (|Stavel ey- Smit h et alTll996|) . In 
brief, an effective bandwidth of 341 MHz, centered around 
1.35 GHz, is divided into 874 frequency channels and sam- 
pled using 2 bits every 64 /is. In April 2009 (MJD = 54939) 
a bright radio pulsar, PSRJl 622-4950, was discovered. It 
has a period (P) of 4.326 s and a dispersion measure (DM) 
of 820 cm~^ pc, corresponding to a distance of ^ 9 kpc, if the 
Cordes-Lazio model for the distri bution of the electrons in 
the interstellar medium is assumed (ICordes & La zio"2002!). 

After the HTRU discovery we embarked on a timing cam- 
paign, observing PSRJ1622-4950 a total of 110 times at 
1.4 GHz and 10 times at 3.1 GHz, obtaining polarimetric in- 
formation on most occasions. The digital filterbank system 
(DFB3) used to create the folded profiles first converts the 
analogue voltages from each polarization channel of the lin- 
ear feeds into digital signals. It then produces 1024 polyphase 
filterbank frequency channels that are folded at the apparent 
topocentric period of the pulsar into 1024 pulsar phase bins, 
and written to disk every 20 s. Four Stokes parameters are 
recorded. To determine the relative gain of the two polariza- 
tion channels and the phase between them, a calibration signal 
is injected at an angle of 45° to the feed probes. The data are 
analyzed off-line using the PSRCHIVE package0 (iHotan et al.l 
2004) and corrected for parallactic angle and the orientation 
of the feed. The position angles are also corrected for Faraday 

^ See |http://psrchive.sourceforge.net| 
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Fig. 1. — Light curve and profile variations for PSRJ1622-4950. The top plot shows the temporal variation of flux density at 1.4 GHz. Archival timing 
data collected within the framework of the Parkes Multibeam Survey ends just before MJD = 54000, and the discovery of PSR J 1622-4950 by the HTRU 
pulsar survey was made on MJD = 54939. The red asterisks show the detections (with errors smaller than the size of the points), and the blue circles indicate 
observations during which PSR J1622-4950 was not detected (down to a limiting flux density of 1.2 mJy). The purple error bar indicates the observation centered 
at 6.3 GHz from the Methanol Multibeam survey that detected PSR J 1622-4950. When converted into a 1.4 GHz flux density, the large error is dominated by 
the uncertainty in the spectral index; the top value of the bar corresponds to a flat radio spectrum and the bottom value is derived from our best estimate of 
the spectral index. The bottom plots show three different pulse profiles from three consecutive observations, taken on 30 June 2009, 16 July 2009 and 23 July 2009. 



rotation through the interstellar medium using the nominal ro- 
tation measure. 

2.2. Archive Mining 

The area of sky containing PSRJ1622-4950 was cov- 
ered previously b y the Parkes Multibeam Pulsar Survey 
(|Manchester et alJ uQQl) but, somewhat surprisingly given its 
large apparent flux density and DM, the pulsar was not de- 
tected in these data. However, that survey did find two other 
radio pulsars in close proximity to PSR Jl 622-4950 and mon- 
itored them for a number of years, namely PSRs Jl 623-4949 
and J 1622-4944 at an angular separation of 1 1' and 7' respec- 
tively. These archival data were reprocessed using the period 
and DM of PSR Jl 622-4950. The first detection was made 
in observations recorded in June 1999 (MJD = 51334), and it 
was then detected in 1 1 observations between October 2000 
(MJD = 51844) and July 2002 (MJD = 52458) (see the top 
panel of Fig. [TJ. However, in the 14 subsequent observations 
preceding August 2006 (MJD = 53975), the pulsar was de- 
tected only once. No further archival data of the two nearby 
pulsars were available after that date. The Methan ol Multi- 
beam pulsar survey at 6.3 GHz (Bates et aHlin prep.|) also cov- 
ered the part of the sky containing PSR Jl 622-4950, and by 
reprocessing the relevant data, the pulsar was detected in an 
observation taken in April 2007 (MJD = 5421 1). 

2.3. ATCA and Chandra Observations 

Observations of the pulsar were made on 8 Dec 2009 
and 27 Feb 2010 with the Australia Telescope Compact Ar- 
ray (ATCA), an east- west synthesis telescope located near 
Narrabri, NSW, which consists of six 22-m antennas on a 6- 
km track. The observations were carried out simultaneously at 
5.5 and 9.0 GHz with a bandwidth of 2 GHz at each frequency 



subdivided into 2048 spectral channels, and full Stokes pa- 
rameters. The source was tracked for 12-hours in each of the 
EW352 and 7 5 OB array configurations. 

Initial data reduction and analysis were carried out with the 
MIRIAD packaged using standard techniques. After flagging 
bad data, the primary calibrator (PKS 1934-638) was used for 
flux density and bandpass calibration and the secondary cali- 
brator (PKS J1613-586) was used to solve for antenna gains, 
phases and polarization leakage terms. 

We also obtained a 20-ks Chandra observation performed 
with ACIS-I on 10 July 2009 to search for an X-ray counter- 
part. Figure[2]shows the resultant radio and X-ray images. Al- 
though there are a number of radio and X-ray sources present 
in the field of view, a highly polarized radio point source is co- 
incident with the brightest X-ray source, CXOU J162244.8- 
495054, which is also the counterpart of a source seen in 
archival ASCA data, AX J162246-4946. We note that CXOU 
J162244. 8-495054 was one of only two X-ray sources in the 
field without an optical/infrared counterpart. We are therefore 
confident that the polarized source seen in the radio is indeed 
the pulsar and that it has an X-ray counterpart. 

3. ANALYSIS, RESULTS AND DISCUSSION 

3.1. Radio Light Curve 

Throughout the HTRU timing campaign, the pulsar was de- 
tected in every observation with a greatly varying flux density, 
which fluctuated by as much as a factor of ~6 within a 24- 
hour period. These flux variations cannot be attributed to in- 
terstellar scintillation and must be intrinsic to the source, as at 
20cm wavelengths, pulsar flux densities ar e reasonably stable 
at such large DMs (IStinebring et al.ll2Q0Ql) . The pulse profile 

^ See |http ://www. atnf .csiro . au/computing/software/miriad | 
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Fig. 2. — Top: Chandra X-ray observations of the region around PSR 
J1622-4950. The source encircled with a 15'' marker, CXOU J162244.8- 
495054, is coincident with the radio point source in the ATCA image. 
Bottom: ATCA radio observations at 5.5 GHz. The encircled source is PSR 
J 1622-4950. South of the magnetar is a faint ring of extended emission that 
appears to be non-thermal and may be a supernova remnant (see text for 
details). 



averaged over each observation (typically 5-10 minutes) also 
changes shape on short timescales, often from day to day. Al- 
though such variability in both the integrated profile and the 
flux density is very uncommon in pulsars, it is a distinctive 
feature of the two magneta rs whose pulsed emi ssion has been 
detec t ed in the radio band dCamilo et ani2007a : Servlak et alJ 
[2QQ7I: iKramer et aL"20Q7'). The light curve for PSR J 1622- 
4950 and three representative pulse profiles are shown in Fig. 

m 

3.2. Timing 

To attain a stable timing solution of a source that is chang- 
ing its pulse profile from observation to observation, special 
care needs to be taken in the timing procedure. In conven- 
tional pulsar timing a standard profile is created and used to 
calculate each pulse's time of arrival (TOA) at the telescope. 
In the case of PSR Jl 622-4950 each pulse profile is different, 
but we noticed that all profiles seem to consist of the same 
components, only differing by their amplitudes (and some- 
times the amplitude for a specific component was indistin- 
guishable from zer o). Software from the PSRCHIVE package 
dHotan et ani2004l) was used to fit gaussians to one of the pro- 
files that showed all the different components and a model 







Fig. 3. — Extract of the P-P diagram, showing the longer period pulsars 
and those with higher inferred magnetic field strengths. Blue dots represent 
ordinary pulsars, red triangles are AXPs, purple squares are SGRs, and 
the black asterisk represents PSR J 1622-4950. The two encircled AXPs 
are the magnetars that are observed to emit radio pulsations. Dashed lines 
correspond to lines of constant characteristic age and dotted lines to constant 
magnetic field strength. PSR J 1622-4950 has a sHghtly higher magnetic 
field strength than the AXPs that are pulsating in the radio band, but clearly 
higher than that of any of the ordinary pulsars. 



was made to describe how the components build up the pro- 
file. This model was then used to calculate TOAs, by letting 
the amplitudes of the different components vary. 

When TOAs for all our observations were calculated, the 
Tempo softwar43 was used to analyze the rotational history 
of PSR Jl 622-4950. The resultant fit requires at least 9 pe- 
riod derivatives to obtain reasonable phase conn ection for all 
our data, which is not uncommon for magnetars ( Camilo et al.l 
|2008|). To monitor variations of the period derivative with 
time, we have also measured P and P over several time in- 
tervals separately. The time spans were chosen to be as long 
as possible while keeping phase connection without adding 
higher period derivatives. This usually resulted in groups of 
TOAs spanning about 40 days each. It appears that P is fluc- 
tuating within a factor of ~2 in time. From this analysis we 
derived an average period derivative, P = 1 .7 x 10~^^ s s~^ The 
parameters from the timing analysis are listed in Table [T] 

The timing solution for PSR Jl 622-4950 implies a very 

high surface magnetic field strength of5 = 3.2xlO^^G \PP 
^ 2.8 X 10^^ G and a ch aracteristic age of r^ = P/(2P) ^ 
4000 years (e.g. lLorimer & Kramer 2005). This is the highest 
surface magnetic field of any radio pulsar known to date. Fig- 
ure[3]shows an extract of the P-P diagram, where the currently 
known magnetars are clustered in the upper right corner, with 
P > 2 s and P>5xlO~^^ss~^ The two radio emitting magne- 
tars are found in the lower parts of the magnetar group, with 
periods that are short in comparison to the other AXPs and 
inferred magnetic field strengths in the lower half of the mag- 
netar range, adjacent to the ordinary radio pulsars with the 
highest surface magnetic fields. The gap between these two 
groups of pulsars is narrowing with new discoveries, such as 
the detection of magnetar-like X-ray bursts from PSR J 1846- 
0258 ( Gavriil et al.. .2008). which was previously thought to 
be solely rotation-powered. PSR J1622-4950 is found in the 



* See http://www.atnf.csiro.au/research/pulsar/tempo/| 
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TABLE 1 
Parameter SUMMARY OF PSRJ1622-4950. 
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Fig. 4.— Polarimetric profiles of PSR J1622-4950 at 3.1 GHz (top) and 
1.4 GHz (bottom) observed at Parkes on the same day (MJD = 55197). The 
blue and red lines represent circular and linear polarization respectively, 
while the black lines represent the total intensity. The full 360° of 
pulse rotation are shown, with the zero point of the pulse longitude set 
at the peak value of the flux density. Absolute position angles at infinite 
frequency are shown after correction for a rotation measure of -1484 rad m~^. 



same region of the plot as the two radio emitting magnetars. 

3.3. Polarization 

The emission from the magnetar is often highly linearly po- 
larized, typically with greater polarization at 3.1 GHz than at 
1.4 GHz (see Fig. [4]). The pulsar also intermittently shows 
significant circular polarization. The Faraday rotation mea- 
sure derived from these observations is very large, RM = - 
1484±lradm~^. A detailed analysis of these polarization 
data will be reported elsewhere but we note the similarity be- 
tween the polarization properties of PSR Jl 622-4950 and the 
other radio pulsating magne tars IE 1547-5408 (Camilo et alj 
l2008i) andXTEJ1810-197 (JKramer et al.ll2007i:iCamilo et all 



3.4. Radio Spectrum 

We combined the various flux density measurements to esti- 
mate the spectral index of the pulsar in the radio band. The av- 
erage flux density at 1.4 GHz is (5'i4oo) = 4.8 mJy (with stan- 
dard deviation (Js^^qq = 2.8 mJy) and at 3.1GHz is (^'3100) = 
4.9 mJy (os^^qq = 2.5 mJy). The ATCA measurements yield 
(^5500) = 13 ± 1 mJy and (^'9000) = 14.3 ± 0.8 mJy, similar to 
the flux density of 12±2mJy obtained from the Methanol 
Multibeam pulsar survey at a centre frequency of 6.3 GHz. 



Parameter 



Valu 



Observed 



Derived 



Right ascension (J2000) ^ 

Declination (J2000) ^ 

Galactic longitude ^ 

Galactic latitude ^ 

Epoch 

Spin period (P) 

Period derivative (P) 

Dispersion measure (DM) 

Flux density at 1400 MHz (^uoo) 

Rotation measure (RM) 



16^22"^44\80(3) 
-49°50'54'^4(5) 
333.85 
-0.10 

MJD 55080 
4.3261(1) s 
1.7(l)x 10-11 ss-^ 
820(30) cm-3 
4.8(3) mJy 
-1484(1) rad m-2 



'pc 



Distance ^ 

Surface magnetic field (B) 
Characteristic age (tc) 
Spin down luminosity (E) 
X-ray luminosity (Lx) ^ 



?^9kpc 

2.8x101^^0 
4kyr 

8.5x10^3 ergs-i 
2.5x10^3 ergs-i 



^The position refers to the coordinates of the X-ray source. 

I' The distance is calculated from the dispersion measure with the Cordes- 

Lazio model. 

'^Assuming the dispersion measure derived distance. 

The spectrum therefore appears to have a positive spectral in- 
dex, highly unusual for pulsars which have a mean spectral 
index of about -1.6 ( Lorimer et al.lll995 ). Again, this peculiar 
spectral behaviour is observed in the other two radio pulsat- 
ing magnetars, which both have a flat (or inverted) radio spec- 
trum (Cami lo et al . 2008 ; L azaridis et a l. 2008; Camilo etaH 
I2007d) . 

3.5. X-ray Emission 

For ordinary rotation-powered pulsars, the X-ray lumi- 
nosity (Lx) is much smaller than the spin-down luminosity, 
E = 471^ I PP~^; on average Lx ^ lO"^^ (iBecker & Trumpet 
Il997|) . In contrast, magnetars are observed to have Lx^E. 
An estimate of the unabsorbed 0.3-lOkeV flux for CXOU 
J162244.8-495054 gives 2.6 x 10"^^ergcm"2 s"^ assuming 
a Galactic hydro gen column density in tha t direction of A^h 
= 2x 10^2 cm-2 (iDickev & Lockmanl[T99Qh and a blackbody 
spectrum with kT = 0.4 keV, typical of a quiescent magnetar. 
This implies Lx(0.3-10keV) ^ 2.5 x 10^^ erg s"^ at the es- 
timated distance of 9kpc and, given L ^ 8.5 x 10^^ ergs~^ 
for this magnetar, Lx ~ 0.3L. Since the distance calculated 
from the DM has a large uncertainty and is known only to 
within a factor of two or so, the derived X-ray luminosity may 
be in error by up to a factor of four. Even within these er- 
rors, the ratio is significantly higher than for ordinary pulsars, 
but well within the range for the magnetars IE 1547-5408 
and XTEJ1810-197 which have Lx ^ O.IL and Lx r^ AE 
(iMereghettil 120081: ICamilo et al.l[2007al: iHalpern et al.][2005h 
in quiescence and larger ratios in outburst. 

The radio emission from XTEJl 8 10-197 was discovered 
immediately following a strong X-ray outburst. It has since 
faded, both in the radio and the X-ray band, and the radio 
pulsations are no longer visible. For IE 1547-5408, the radio 
emission is also highly variable and appears to be 'revived' in 
the periods after its X-ray outbursts. PSR Jl 622-4950 on the 
other hand, has had at least two episodes of non-detections in 
the radio band lasting hundreds of days followed by periods 
of bright radio emission (see Fig. [T]). In addition to the new 
Chandra observation, we searched archival data from Chan- 
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dra, XMM Newton, Rosat, ASCA, Beppo-SAX, Rossi-XTE and 
Swift for an outburst, however no evidence for X-ray flux vari- 
ability and no X-ray outburst at the level of the outbursts seen 
in XTE 1810-197 and IE 1547-5408 in connection to the ra- 
dio pulsations (> lO^^ergs"^) were found since at least as 
early as 2005. It is possible therefore that an enhancement 
of X-ray activity is not a requirement for pulsed radio emis- 
sion by magnetars, however, given the duty cycle of sensitive 
X-ray observations of the field containing PSRJl 622-4950, 
we cannot constrain the occurrence of fainter X-ray enhance- 
ments of the source. What is instead certain is that the ob- 
served X-ray emission from PSRJl 622-4950 is at variance 
with what is observed for the other two radio pulsating mag- 
netars. 

3.6. Possible SNR Association 

If the true age of PSR Jl 622-4950 is similar to its character- 
istic age of 4 kyr, we might expect to see a supernova remnant 
(SNR) surrounding the pulsar. Indeed, 5 of the 9 AXPs and 
at least 1 of the 5 SGRs are located within SNRs (Mereghettl 
120081: iGaensler et al. 2001). Inspecting the ATCA image in 
Fig. O we see a ring of emission centered ~2' south of the 
pulsar location. This ring lacks an infra-red counterpart and 
appears to be non thermal, whereas the extended radio source 
to the south of the ring is clearly thermal in nature. Could 
the ring be the SNR and the pulsar has escaped its bounds? 
If we assume a distance of ~9 kpc to the magnetar and fur- 
ther assume it was born in the centre of the ring, the magnetar 
would need a velocity of ~ 1 300 km s~^ to reach its current lo- 
cation whereas the ring itself would have a lower expansion 
velocity. Such a velocity is high (though not impossible) for 
pulsars but rather low for expanding SNRs. Although the link 
between the ring and the magnetar is a possibility we consider 
it unlikely. 

4. CONCLUSIONS 

The HTRU survey has discovered a radio-luminous pulsar, 
which is highly polarized, has an inverted spectrum, and is 
highly variable in both its pulse profile and flux density. The 



radio pulsar has a faint X-ray counterpart that appears to be 
stable in flux, with a value that is typical of a quiescent magne- 
tar. The pulsar shares many of the properties of the two known 
radio magnetars and we therefore conclude that PSR J 1622- 
4950 is indeed a magnetar, the first discovered through its ra- 
dio emission. This discovery not only adds a new member to 
the magnetar family, but also highlights unprecedented fea- 
tures of the emission of the magnetars across the electromag- 
netic band. At odds with what is observed in other sources, 
PSR Jl 622-4950 indicates that bright radio emission can be 
present even when a magnetar displays an X-ray luminos- 
ity typical of a quiescent state. Moreover, PSR Jl 622-4950 
shows that pulsed radio emission can either exist without the 
occurrence of a strong X-ray outburst, or occur a long time 
(> 5 years) after the outburst. Alternatively, the radio pulsa- 
tions could be triggered by a modest increment of X-ray ac- 
tivity, that escaped detection in this case. We finally note that 
the extreme variability in the flux density of PSR J 1622-4950 
also demonstrates the advantages of surveying the radio sky 
at regular intervals with even modest sensitivity. This high- 
lights the potential of the upcoming radio facilities like the 
LOR\R, ASKAP or the SKA which promise to characterize 
the dynamic radio sky at an unprecedented level. 
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